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Dissipative effects during neutrino decoupiing in the eariy universe create a smali backreaction on 
the Hubble rate, and lead to a small rise in temperature and entropy. We use a simplified thermo- 
hydrodynamic model, which provides a causal approximation to kinetic theory, in order to estimate 
the backreaction effects and the entropy production. 
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I. INTRODUCTION 



Non-equilibrium processes in the early universe are typically associated with dynamical transitions or particle de- 
couplings. In the case of neutrino decoupling, the standard approach is to treat the process as adiabatic (see e.g. [jf)). 
The small non-equilibrium effects are thus usually neglected, which provides a reasonable approximation. However, 
given the increasing accuracy of cosmological observations and theoretical modeling, it is worthwhile revisiting the 
£S| ■ standard equilibrium models of processes such as neutrino decoupling, in order to see whether non-equilibrium cor- 
^ ' rections can lead to observable consequences. Recently, non-equilibrium corrections in neutrino decoupling have been 
calculated in a number of papers, using complicated kinetic theory and numerical computations (see Q for a short 
review). The corrections are very small, as expected. For example, in |^-|5|] it was found that non-equilibrium effects 
lead to a small change in the decoupling temperature for neutrinos. Spectral distortions have also been analyzed ||, 
showing the remarkable fact that they amount to as much as 1% or more for the higher-energy side of the spectrum. 
Although these corrections in the spectrum, energy density and temperature of the neutrino component have hardly 
^\ . any effect on primordial helium synthesis, yielding a change in the mass fraction of ~ 10 -4 , they can lead to other 
effects that may be observable. Thus it is shown that the non-equilibrium increase in neutrino temperature, which 
leads to an extra injection of energy into the photon spectrum, leads to a shift of equilibrium epoch between matter 
and radiation which, in turn, modifies the angular spectrum of fluctuations of the cosmic microwave background 
radiation @,|). 

Despite the accuracy of these models in obtaining corrections to the decoupling temperature and distribution 
function due to non-equilibrium effects, they still make use of the standard Friedman equations for a perfect (i.e non- 
dissipative) fluid. This leads to the physically inconsistent situation in which, say, the energy density and expansion 
evolve in time like a radiative fluid in equilibrium. One expects that small distortions in the particle equilibrium 
distribution function should be reflected in the macroscopic (i.e fluid) description, as given by the stress-energy tensor, 
by adding a bulk viscous pressure to the equilibrium one. Here we consider an alternative thermo-hydrodynamic model 
of dissipative effects in neutrino decoupling, simple enough to produce analytic solutions for the backreaction effects 
on the universal scale factor, and estimates for the entropy production due to dissipation. As explained above these 
effects are not the focus of recent papers, which use sophisticated kinetic theory models focusing on the neutrino 
temperature. Our simplified approach cannot compete with these models for accuracy and completeness, but it has 
the advantage of simplicity, allowing for a qualitative understanding of effects not previously investigated in detail. 
A similar approach has previously been developed in || to the reheating era that follows inflation. 

The thermo-hydrodynamic model is based on an approximation to kinetic theory which respects relativistic causality. 
This approximation is the Grad moment method, leading to the causal thermodynamics of Israel and Stewart |lC| ] 
in the hydrodynamic regime (see also ]i"l| ] for an alternative but equivalent approach). This causal theory is a 
generalization of the more commonly used relativistic Navier-Stokes- Fourier theory. The latter, due to Eckart jL2j, 
may be derived via the Chapman-Enskog approximation in kinetic theory. The resulting theory is quasi-stationary 
and noncausal, and suffers from the pathologies of infinite wavefront speeds and instability of all equilibrium states 
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|L3|]. The main new ingredient in the causal transport equations is a transient term which contains the relaxation 
time. Our simple model is based on a one-component fluid. In ||l4| , relaxation time processes are incorporated in 
a two-fluid model. In this setting, electrons and positrons on the one side and neutrinos and antineutrinos on the 
other side, are found to be in two different equilibrium states with slightly different temperatures. The system evolves 
towards a state of thermal equilibrium in a characteristic relaxation time. 

Dissipative effects in the decoupling of a given species of particles arise from the growing mean free path of the 
decoupling particles in their weakening interaction with the cosmic fluid. Eventually the mean collision time exceeds 
the gravitational expansion time, and decoupling is complete. A hydrodynamic model may be used to cover the early 
stages of the decoupling process, but it will eventually break down when the mean collision time becomes large enough 
S 

In the conditions prevailing at the time of neutrino decoupling, it is reasonable to neglect sub-horizon metric 
fluctuations and treat the spacetime as a Friedmann model. (The incorporation of perturbations in our model would 
use the covariant formalism for dissipative fluids developed in |l6j].) The dynamical effects of spatial curvature and any 
surviving vacuum energy will be negligible, so that we can reasonably assume a spatially flat geometry. Furthermore, 
we assume that the average 4- velocities of the neutrinos (regarded as massless) and of the photon-electron-positron gas 
are the same. With all these assumptions, only scalar dissipation is possible. Dissipation during neutrino decoupling 
arises because the falling temperature lowers the interaction rate with leptons as the lepton mass can no longer be 
ignored relative to the thermal energy. Thus dissipation is directly reflected in a deviation of the equation of state 
from the thermalized radiation form p = ^p. Within a hydrodynamic one-fluid model, such dissipation is described via 
bulk viscosity, which vanishes in the p = ^p limit, but is nonzero otherwise. We will use the full (i.e. non-truncated) 
version of the causal transport equation for bulk stress. 



II. CAUSAL TRANSPORT EQUATION FOR BULK STRESS 

The particle number 4-current and the energy-momentum tensor are 

iV a = nu a , T ab = pu a u b + {p + U)h ab , 

where p is the energy density, p is the equilibrium (hydrostatic) pressure, n is the particle number density, II is the 
bulk viscous pressure, and h ab — g ab + u a u b is the projector into the comoving instantaneous rest space. Particle and 
energy-momentum conservation 

\J a N a = , V b T ab = , 

lead to the equations 

h + 3Hn = Q, (1) 
p + 3H(p + p + U) =0, (2) 

where H is the Hubble expansion rate. The specific entropy s and the temperature T are related via the Gibbs 
equation 

nTds = dp- ^^-dn . (3) 
n 

Then it follows that 

nTs = -3HU , (4) 

where II is always non-positive. The Grad moment approximation in kinetic theory (or phenomenological arguments) 
leads to the full causal transport equation JlGj] for If: 



7-ri + n = -3(H - irn 



t C T 



(5) 



where r is the relaxation time scale, which allows for causal propagation of viscous signals, and C < is the bulk 
viscous coefficient as given below. Quasi-stationary, noncausal theories have r = 0, which reduces the evolution 
equation (||) to an algebraic equation II = —3(H. This leads to instantaneous propagation of viscous signals. Note 
also that the causal relaxational effects lead to a small increase in the sound speed over its adiabatic value [jl7f : 
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(6) 



This result, which is not well known, is derived in the appendix. 

The approximation used in deriving the transport equation (also in the quasi-stationary case) requires that |IT| <C /?, 
which is reasonable for most dissipative processes (see p8[ for a nonlinear generalization of the causal transport 
equation.) 

Equation (|^) as it stands is known as the full or non-truncated transport equation for bulk viscous pressure [^| ^| . 
When the term containing the square bracket on the right is neglected, we get the truncated equation which is usually 
used. Under many conditions, truncation leads to a reasonable approximation. We will use the full equation. 

Taking n and p as independent variables, the Gibbs equation (0) leads to the integrability condition 



dT 



(7) 



and together with the energy conservation equation 



this gives the temperature evolution equation 
SHU. 



(8) 



The first term on the right accounts for adiabatic cooling due to expansion, whereas in the second term, viscosity 
contributes to heating of the fluid (note that II is always non-positive). 
Using equations ([|) and (||), the Gibbs equation takes the form 



n 2 Tds 



n5HU 



3H(p + p) + 3HU 



,1,,+ (/- f/M ( ^ 



— dp — dp 
P 



(9) 



As expected we learn from the last equation that when the fluid is perfect (II = 0), the specific entropy is conserved 
along the flow lines (s = 0). Furthermore, if a barotropic equation of state for n holds, i.e. n = n(p), then ds = so 
that s is a universal constant, the same on all flow-lines, and the fluid is called isentropic.[] Yet, as Eq. (|^) shows, 
this is no longer true in the presence of dissipation, i.e. a barotropic particle number density no longer forces ds to 
vanish. 

For simplicity, we assume the linear barotropic equation of state 



(10) 



where 7 is constant and we are interested in the case 7 



The adiabatic speed of sound c s is given by 



dp 

dp 



which for a perfect fluid (either barotropic or not) becomes 



When Eq. ( |lO| ) holds then c s = ^7 — 1. Using Eq. (|TTJ) and the integrability condition (j^), we find 



j 1 = p(i~ 1 )/i p 



nV7 



(11) 



where F is an arbitrary function which satisfies F — 0. If T is barotropic, then F is constant and we have a power-law 
form with fixed exponent for the temperature |l7| , |22| 



T oc p 



( 7 -l)/7 



(12) 



1 The same reasoning applies when the temperature is barotropic. 
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In the non-dissipative case, these barotropic equations for p and T are compatible with the ideal gas law 

p = nT, (13) 

but in the presence of dissipation this is no longer true. In effect, equations (|l(i|), ( |l~2| ) and ( [l3| ) imply n cx p 1 ^, i.e. 

h lp 
n j p 

which implies, by using Eq. (||), that II = 0. We shall drop in the sequel a barotropic equation of state for the 
temperature in favour of the more physically appealing equation of state ( |l3| ) together the 7-law in ( |l0| ) . 



III. DISSIPATION IN NEUTRINO DECOUPLING 



A hydrodynamic approach in the expanding universe requires a particle collision time t c short enough to adjust to 
the falling temperature. As the natural time-scale for the expanding universe is H , we have 

t c < H- 1 . 

If t c <C i/ -1 , then an equilibrium state can in principle be attained. Dissipative phenomena could play a prominent 
role for t c ~ H~ l . 

We learn from kinetic theory that t c is determined by 

t c = — , (14) 
nav 

where n is the number density of the target particles with which the given species is interacting, a the cross-section 
and v the mean relative speed of interacting particles. For the decoupling of massless neutrinos in the early universe, 
v = 1, the target number density is that of electrons, and p3| ] 

where G F is the Fermi coupling constant. At the neutrino decoupling temperature Td, we have m /Td w |, so that 
the rest mass energy m c of electrons starts to become important. Since the electron number density in the radiation 
dominated era evolves as n c cx a~ 3 , where a is the scale factor, we have from Eq. (|lj) that 

*c«^- (15) 

Dissipation due to massless particles with long mean free path in a hydrodynamic fluid is described by the radiative 
transfer model. The bulk viscous coefficient takes the form [Q 

C = 4rT 4 r 2 i c , (16) 

where r is g times the radiation constant and T measures the deviation of p/p from its pure-radiation value: 

r =Ht), 

where p and p refer to the pressure and energy density of the radiation/matter mixture as a whole. Since we assume 
the linear equation of state (|l0|), it follows that T is a perturbative constant parameter in our simple model: 

r = |- 7 «i. 

The assumption that T is constant relies on the assumption that decoupling takes place rapidly. Since standard 
adiabatic treatments of decoupling [jjj assume instantaneous decoupling, this assumption should be a reasonable first 
approximation. 

We may neglect the — 3(H term on the right of the transport equation (||), since it is 0(T 2 ). Note that our simple 
model would thus break down in the quasi-stationary Eckart theory, since it would immediately lead to II = 0(T 2 ). 
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The relaxation timescale r in causal radiative transfer [B5| is given by r = t c . The term ('/('on the right of Eq. 
becomes 



= h + o(r) : 



on using equations (||) and ([l5]). The full transport equation (||) becomes, to lowest order 

rri + n = -4ri?n . 



(18) 



(We can think of the right hand side as an effective source term relative to the truncated transport equation.) We 
can rewrite this in the standard truncated form as 



T*n + IT = , 

where the effective relaxation time acquires an expansion correction: 



1 + AtH 



(19) 



(20) 



The amount of reduction depends on the size of r = t c relative to H. The hydrodynamical description requires 
tH < 1. If tH <C 1, then r* w r. But if tH is close to 1, the reduction could be significant. 
The Friedmann equation 



p = 3H 2 , 

together with Eq. (§) leads to 

n = -2H- (4-3T)H 2 . 
On using equation (^) we get from (|l8| ) the evolution equation for H 



where 



H + HH(8 - 3T + N) + H 3 (2 - |r)(iV + 4) = 0, 



N=(tH)-\ 



(21) 
(22) 
(23) 
(24) 



which is of the order of the number of interactions in an expansion time. Now, from equations (|l0|), (p^), (15) and 
(pi|) we have 



N = 



Ha 
HdddJ 



(25) 



where the expression n oc a 3 has been used and ad and = -ff(ad) are the values at which N = 1, so that is 
determined by the equation 



t c (a d )H(a d ) = 1 , 



(26) 



Changing the independent variable to the scale factor a, developing equation (|23|) and collecting the previous results, 
yields 



a 2 HH" + a 2 H' 2 + aHH' 



(2 - f T) H 2 



9-3r 



Ha 



f— Y 

\H d a d J 



= 0, 



where a prime denotes d/ da. We expand H as 

H = H + SH where 5H = Th + 0(T 2 ) . 



(27) 



(28) 
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The equilibrium Hubble rate H corresponds to the thermalized radiation state p — so that T — 0, and Eq. (51 
becomes 



a 2 HH" + a 2 H' 2 



9aHH' + 8H 2 + (aHH' + 2H 2 ) ( -^-] = . 

\H d a d J 



The unique power-law solution is the well-known perfect radiative solution 



«o 
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where ao marks the start of the dissipative decoupling process, so that H = H for a < a^. 
Substituting Eq. (pq) into ( |2^ ) and using the fact that 

H d a d a Q 

we find that to 0(T): 



(29) 



2, II 



a h 



(?) 



h! 



4 + 2^ 



v a 



h = ^Hq 



0£ \ /Od 

a d J \ a 



Defining a = a/a d , we can rewrite this as 

d 2 h 
da 2 



"5 


1 " 




' 4 


2 " 


- - 

a 




da 


a 2 


h — 



h = 



2 H O a o) tf- 



Now we use the following general result [£6| : if ip is a solution of 

V" + f(x)y' + g(x)y = k(x) 

when k — 0, then the general solution is 



y = Ci<p + C 2 f 



dx 



tp 2 E 



ipEkdx I dx , 



(30) 



(31) 



where E — exp J fdx. By inspection, a solution of the homogeneous equation (^Tj) is I /a 2 . It follows that the general 
solution is 



(32) 



where ci and C2 are arbitrary integration constants and Ei is the exponential-integral function 



— dv = C + In a; + ■ 



a; 



with C denoting Euler's constant. 

By equations (0) and (32), the bulk stress to first order is 

II = (3H 2 - 4Hh - 2h'Ha)T, 



(33) 



This expression holds for a > ao, where ao marks the onset of dissipative evolution. Thereafter, the bulk stress 
decays according to the causal law (|l9|). In order to relate the constants c\ and C2, we require, according to standard 
matching conditions, that H is continuous. Thus h(ao) = 0, which fixes c\\ 



ci = -c 2 Ei 



1 / Od 

3 Uo 



In 



(34) 



Thus, using Eq. (p2|), we see that the backreaction of the dissipative decoupling process on the expansion of the 
universe is given by 



G 



1 /Od 



3 Vet 



Ei 



3 V a 



SH = H jc 2 ^Ei 

+ o(r 2 ). 

Substituting Eq. (]3~I|) into Eq. ((33|), we find that the bulk stress becomes 



2 \a 



II = p < 2c2 exp 



1 (OdV 



3 V a / 



r + o(r 2 ), 



(35) 



(36) 



where p = 3H 2 is the equilibrium energy density. Since II < 0, we require C2 < 0. Below we find a prescription for Ci 
in terms of physical parameters. 



IV. CONCLUSION 



In order to complete the model, we need to determine the remaining arbitrary constant C2 in terms of physical 
parameters. A rough estimate, which is consistent with the simplicity of the model, arises as follows. We estimate 
the duration of the dissipative process as 



Aa 



fld — do , 



(37) 



i.e. we assume that the process ends at a<j- Then by Eqs. (| 
to decoupling is approximately 



and (fL3|) , the fractional viscous rise in temperature due 



AT 
~T~ 



II(oo) Aa 

p(a ) a 



(38) 



We can consider the fractional temperature increase as an input from previous kinetic-theory investigations (as 
described in the introduction), which typically predict it to be O(1O _3 ).0 Then equations (|36|)-([38"|) and ( |3^ ) allow us 
to estimate the constant ci in terms of the physical parameters a^/ao, AT /T and T as: 



c 2 r 



1 AT 


exp 


_I ( 2d\ 

3 \a J 


2~ i 


1 1 





(39) 



Finally, we can also estimate the entropy production due to decoupling. By Eqs. 
in entropy per particle is approximately 



A> 



AT 



and (|38|), the viscous increase 



(40) 



Our model describes the response of the cosmic fluid to a bulk stress, which is a very simple thermo-hydrodynamic 
approximation to more realistic kinetic theory models of neutrino decoupling, but which nevertheless accommodates 
the dissipative effects and respects relativistic causality. The simplicity of our model allows us to derive analytic forms 
for the dynamical quantities and the backreaction effects, but it does not incorporate a mechanism for bringing the 
dissipative process to an end. 
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APPENDIX A: CHARACTERISTIC VELOCITIES FOR BULK VISCOUS PERTURBATIONS 



Following [fnl], we derive equation (^|) for the dissipative contribution to the sound speed. The full analysis of the 
causality and stability of the Israel-Stewart theory was performed in a series of papers by Hiscock and Salmonson 
p3| , p8| . They showed that both issues are closely related and obtained general expressions for the characteristic 
velocities for dissipative perturbations. Here we extract from their general expressions specific results for the case in 
which only bulk viscosity is present. 

The purely bulk viscous case stems from the general expressions of by setting all the coefficients coupled to 
heat flux and shear viscosity to zero. This yields for the speed of propagating transverse modes 



(p + p)a 2 1 +2a 1 + (3 1 
2fo\p 1 {p+p)-l] 



0, 



which is what one expects for scalar sound-wave perturbations. Equation (128) of |13] ] governing the speed v = v L of 
propagating longitudinal modes becomes, on dividing by /3o/?2 and setting ao — a\ — 0, 



~2n 



T 
1 

+ rtf 2 

Dividing by (3\ and taking /3\ — 



dT\ 
dn) . 



(P + P) 



nT 2 
(P + P) 
00, we have 



dT 
1h 



Pi\(p + P) 



dp 
dp 



1 

Jo 



dT 



dp 

dp 



dp 

dp 



1 

A) 



dT 
dn 



1 



(Al) 



The first term on the right is the adiabatic contribution c 2 to v 2 , and the second term is the dissipative contribution 
c^, which, requiring v 2 < 1, leads to 



cl- 



(p + p)r 

We also learn from |13) that causality and stability require 



< 1 



n 3 (\) = (p + P )\i-\ 



dp 
dp 



{p+p)Po 



< 0. 



(A2) 



(A3) 



for all A such that < A < 1. This condition is seen to be hold on account of the inequality ( A2). 

The expression for Cb refines and corrects the statement in J2i| (the first paper to apply causal bulk viscosity in 
cosmology) that C/ P T — 1 1S required by causality. 
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